Ortho-para transition in molecular hydrogen 
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The radiative ortho-para transition in the molecular hydrogen is studied. This highly forbidden 
transition is very sensitive to relativistic and subtle nonadiabatic effects. Our result for the transition 
rate in the ground vibrational level r(J = 1 — > J = 0) = 6.20(62) • 10"^* year~^ is significantly 
lower in comparison to all the previous approximate calculations. Experimental detection of such 
a weak line by observation of, for example, the cold interstellar molecular hydrogen is at present 
unlikely. 
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The hydrogen molecule in the ground electronic state 
can exist in a nuclear triplet state {S = 1, ori/io-H2) 
with the odd angular momentum L, or in a singlet state 
{S = 0, para-H2) with the even L. The question we rise 
is, what are the physical mechanisms for possible transi- 
tions between these two classes of states. The nonradia- 
tivc transition, for example, in the interstellar molecu- 
lar hydrogen is mostly induced by collisions with atomic 
H. The corresponding rates were obtained by Sun and 
Dalgarno in The radiative transition which is much 
weaker can in principle take place at sufficiently low den- 
sities and temperatures. The relativistic spin-orbit in- 
teraction (nuclear spin and the electron momenta) is the 
most obvious source of this transition as it mixes slightly 
the 0-H2 and P-II2 states, see Eq. (jlj. This effect has 
been considered in the original work of Raich and Good 
in 0, although not in a complete and systematic way. 
It has happened that a tiny nonadiabatic correction to 
the total H2 wave function significantly changes the theo- 
retical predictions for this rate. Moreover, the spin-orbit 
mixing is not the only one effect, which makes this transi- 
tion possible. There are also relativistic corrections to the 
El coupling to the electromagnetic field which is barely 
known in literature. These corrections in the context of 
H2 molecule have been derived for the first time by Do- 
delson in Q using the Feinberg-Sucher formalism. Here 
we rederive this result in a much simpler way. Because of 
the summation over the infinite H2 spectrum, the calcu- 
lations of ortho-para transition amplitude are not com- 
pletely trivial. The most elaborate so far calculations by 
Raich and Good including Dodelson corrections Q 
gave the rate of 1.85(46) • 10"^'^ year~^, which did not 
include all important contributions. The purpose of this 
work is to present a complete theoretical description of 
the radiative ortho-para conversion of H2 molecule, in- 
cluding the results of direct numerical calculations of the 
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transition probability for the lowest rotational levels. 

The interaction of an arbitrary molecule with the 
electromagnetic field, whose characteristic wavelength is 
much larger than the size of this molecule, including the 
relevant spin-orbit interaction is [3| 
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In the above, A, b are the indices of the nuclei and 
electrons respectively, XA,Xb are the coordinates of the 
nuclei and electrons with respect to the mass center, 
XAb = XA — Xb, and gA is the nuclear g-factor. Moreover, 
the electromagnetic fields E, B and their derivatives are 
assumed in Eq. ^ to be at the mass center. 

For 0-H2 (the first excited rotational state) the nuclear 
spin 5=1 couples to the orbital angular momentum 
J = 1 of the nuclei, giving the total angular momentum 
characterized by quantum numbers F = 0,1,2. In the 
P-II2 (the ground rotational state) the total angular mo- 
mentum is F = 0, therefore the one photon ortho-para 
transition from any other F = level is strictly forbid- 
den, while F = 1 level decays by El transition and F ~ 2 
level decays by the M2 transition. 

Let us first consider the M2 transition from F = 2 level 
of 0-H2, to F = level of P-H2. This transition comes 
from the following interaction with the electromagnetic 
field, which is obtained from Eq. ^ 
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with TZ = XA — XB and the proton g-factor gp — 
5.585 694 713(46) [5]. From this Hamiltonian one obtains 
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the transition rate T{M2) = r2 



r2 = 2a^^J2 
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where TZq is the average distance between protons, and 
the last equation holds for the nuclear H2 wave function 
which is strongly peaked around TZq, as it is the case for 
the lowest rovibrational levels. 

The calculation of the rate for El transition from 
F — 1 level of 0-H2 to F = level of P-H2 is more 
complicated as it includes corrections to the wave func- 
tion coming from the spin-orbit interaction. Since it 
is AF — 1 transition, the operators in the interac- 
tion Hamiltonian in Eq. ([T]) can be simplified, namely 
e*-''^ {xA y-SA)^ 12 and this Hamiltonian becomes 
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The resulting transition rate r(i?l) = Fi is 
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where 
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and H is the 4-body nonrelativistic Hamiltonian. In or- 
der to simplify the evaluation of Q'^ we perform the adi- 
abatic expansion, namely the expansion of the resolvent 



in the kinetic energy of the nuclei 
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similarly for ^^_b -i ^'^'^ the expansion of the wave func- 
tion 

0E = 7/^(flc,:E2c;7^)Ao(7^)/\/4^ + <50s, (12) 
4>'h = V^(flc,2?2c;7^)A^(7^)/V4^+(5,/)^, (13) 

with = Ai T&ln and with normalization 

j dTzn^ xliiz) = j dnn^ xHn) = i . (i4) 

While the exact nonadiabatic correction is unknown, we 
need only the first order m/mp part of the correction, 
which explicitly depends on the nuclear state Q 
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We introduce now the perturbed electronic wave func- 
tions 



1 


Ebo 


— Hbo 




1 


Ebo 


— ^BO 




1 


Ebo 


-i^BO 



(xlc+4c)V', (17a) 
(17b) 

V^V, (17c) 



to simplify the matrix elements of in Eq. (jTU]) 
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with n = TZq/TZq. One notes that derivatives of nonhnear 
and Unear parameters in ij;, see Eq. (|24l) . with respect to 
TZ do not contribute to the above matrix elements, which 
significantly simplifies the numerical computations. 

Results for Xi can be expressed in terms of dimension- 
less factors Fi 
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which are chosen in such a way, that Fi (m a TZ) vanish 
at 7?. = and approach 1 for 7?. — > cxd. However, in the 
case of F2 this large TZ limit is only a rough approxima- 
tion, since we have not been able to perform this limit 
analytically. We can now return to Eq. (|18p and obtain a 
compact formula for the matrix element in the transition 
rate Ei of Eq. jS]) 
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Numerical evaluation of the E2 rate according to for- 
mula ([3|) is straightforward. To obtain the ortho-para 
energy spacing to and the average internuclear distance 
TZo , we employed the accurate Kolos-LeRoy-Schwartz in- 
teraction potential which includes the adiabatic and 
relativistic energy corrections. With this potential we 
solved numerically the radial Schrodinger equation to ob- 
tain the energies and wave functions corresponding to the 
lowest ortho and para levels. The numerical values used 
here are uj — 27t ■ 118.49 cm^^ and TZq = 1.449 au, and 
the resulting M2 transition rate with physical constants 
from Ref. 5] is 



T2 = 1.07(1) • 10"^"* year 
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The accurate evaluation of Fi and the corresponding 
electronic matrix elements Fi in Eqs. (pij) is a challeng- 
ing task. We have represented the electronic ground state 
wave function as well as the first order perturbed func- 
tions defined by Eqs. p7p . in the form of properly sym- 
metrized linear combinations ^ = CkPg.u4'k of Gaus- 
sian geminals 

(pk = Sfe (24) 
X exp (-Qffe xjj^ - f3k xlg - Ck - Vk xlg - •jk . 

The projection operators 



Pg,u = 7(l+^i2)(l±*) (25) 



ensure the proper symmetry with respect to the exchange 
of the electrons and with respect to the inversion oper- 
ation, yielding singlet gerade or ungerade functions. Re- 
quired S , S , or n symmetry of the electronic wave 



function was imposed by the Cartesian prefactor . The 
linear and the nonlinear parameters were optimized vari- 
ationally with the goal function being the ground state 
energy in the case of the unperturbed wave function ijj or 
pertinent Hylleraas functional 
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in the case of the perturbed functions V'fe- Table [J shows 
explicitly the elements defining particular functions with 
the assumption that the molecule is placed along the 
Cartesian X axis. The unperturbed wave function has 

TABLE I: The definitions of the functions used in the com- 
putations. 
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been expanded in 600-term basis set which enables the 
electronic ground state energy to be obtained with an 
error of only 3 • 10~^ a.u. A 1200-term expansions have 
been employed to represent the perturbed functions. Val- 
ues of the J" functionals corresponding to the optimum 
parameters are displayed in Table |TT1 

The general formulas ([20]) . in the particular case of 
the molecule oriented along the X axis, can be explicitly 
written as follows 
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Table ITIl contains all the expectation values appearing in 
Eqs. ((27l) as well as the final Xi and Fi values computed 
at TZq — 1.449 bohr. To check the correctness of our 
codes we performed additional calculations at large in- 
ternuclear distance [TZ = 12.0 bohr) and compared the 
resulting goal functions and the expectation values with 
analytically derived asymptotic values. This comparison 
is presented in Table |TT1 

Using Eqs. Q, ^F^, and the Table [ll] one obtains the 
numerical value for the El transition rate 



Fi = 1.68(17) • 10"^^ ycar-i , 
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and finally the rate averaged over the total angular mo- 
mentum F 



F = (5E2 + 3Ei)/9 = 6.20(62) • 10"^'' year 



(29) 
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TABLE II: Numerical values of the optimum goal functions 
and the expectation values comprising the Xi factors. The 
asymptotic value for F2 is approximate. 

7^(, = 1.449 7^ = 12.0 Asymp. 



Our result for the averaged transition rate is in dis- 
agreement with the result of Dodelson [a], r = 1.85(46)- 



10"-'^^ year" ""^j which is in turn based on the former work 
of Raich and Good [2,] and included direct coupling of nu- 
clear spin to the radiation field. We confirm in this work 
the existence of these additional couplings, which here 
are expressed by the 3rd term in Eq. ([T]). In our opinion, 
the difference from our result is due to the omission of 
the M2 transition, the omission of the nonadiabatic con- 
tributions corresponding to X2 and in Eq. (|20p . less 
accurate w, and due to a lower accuracy of the numerical 
calculation of the matrix elements in Ref. In partic- 
ular, without X3 the overall rate F would be about 24% 
larger. 

The possibility of the experimental detection of the 0- 
p H2 line is questionable. Much stronger E2 lines have 
already been observed at the Infrared Space Observatory 
ISO and served for estimation of the temperature of in- 
terstellar hydrogen clouds and of the ratio of abundance 
0-H2 to P-H2, which sometimes differs significantly from 
the equilibrium one 0]. The much weaker E\ line has 
not been observed yet. In fact there is a potential op- 
portunity related with the Herschel Space Observatory 
to be launched in 2008 [9!]. Its spectral range covers the 
o-p line at 84.4 /xm, but its resolution is, probably, not 
high enough at this wavelength. 
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